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Surface Enhanced Raman Spectroscopy (SERS) is a well-established technique for en¬ 
hancing Raman signals.®^^ Recently photonic integrated circuits have been used, as an 
alternative to microscopy based excitation and collection, to probe SERS signals from ex¬ 
ternal metallic nanoparticles.^i^HIIl However, in order to develop quantitative on-chip SERS 
sensors, integration of dedicated nanoplasmonic antennas and waveguidesli^^ is desirable. 
Here we bridge this gap by demonstrating for the first time the generation of SERS signals 
from integrated bowtie nanoantennas, excited and collected by a single mode waveguide, 
and rigorously quantify the enhancement process. The guided Raman power generated by 
a 4-Nitrothiophenol coated bowtie antenna shows an 8 x 10^ enhancement compared to the 
free-space Raman scattering. An excellent correspondence is obtained between the theoreti¬ 
cally predicted and observed absolute Raman power. This work paves the way towards fully 
integrated lab-on-a-chip systems where the single mode SERS-probe can be combined with 
other photonic, fluidic or biological functionalities. 

A schematic of the device under study is shown in Figure [TJ a). The fundamental TE-mode of a 
silicon nitride (SiN) rib waveguide excites a periodic array of gold bowtie antennas coated with a 
4-Nitrothiophenol (NTP) monolayer. The pump wavelength for all experiments is set to Ap = 785 
nm and NTP Stokes light (at Xs) is subsequently collected back into the same waveguide mode. 
Fabrication details can be found in the Methods section and a description of the measurement 
setup is outlined in the Supplementary Information SI. A scanning electron microscope image of 
the functionalized waveguide, with cross-sectional area of 220 x 700 nm^, is depicted in Figure 
[^b). Raman spectra of an uncoated and coated waveguide functionalized with 40 antennas are 
shown in Figure [^c). The spectral regions where an NTP Stokes peak is expected (1,080, 1,110, 
1,340 and 1,575 cm^^)l^are highlighted by the cyan shaded areas. Before coating no NTP peaks 
can be distinguished from the inherent SiN background. The peaks at 1,250 and 1,518 cm ^ 
(marked by the black dashed lines) are attributed to interference effects of the Au array which 
act on the scattered background light, and they are also observed on the extinction curves of the 
functionalized waveguides (see Supporting Information S2). Hence they do not represent specific 
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Raman lines. After eoating, four additional peaks appear and eoineide with the expeeted NTP 
Stokes peaks. This demonstrates that SERS signals from single monolayer eoated antennas ean be 
effieiently exeited and eolleeted by the same fundamental waveguide mode. 

Subsequently the dependenee of the SERS signal on the position of the plasmon resonanee was 
investigated to verify that it ean be attributed to a resonanee effeet and not to eoineidental surfaee 
roughness. To this end, waveguides funetionalized with a fixed number of antennas but varying 
bowtie geometries were eonsidered. The relevant bowtie parameters are its length L, gap A and 
apex angle a (Eigurej^a)). By ehanging the length, the antenna resonanee ean be tuned (Li = 90 
nm, L 2 = 115 nm and L 3 = 140 nm for fixed a = 60° and A = 40 nm). Extinetion speetra are 
plotted in Eigure|^b) while the eorresponding Raman speetra are depleted in Eigurej^e). The 
Raman speetrum of a referenee waveguide without any Au funetionalization is also shown. Even 
after eoating the referenee waveguide does not generate NTP peaks, so any Raman signal indeed 
originates from the antenna region and does not eontain eontributions from spontaneous Raman 
seattering along the waveguide.^ The Li resonanee is detuned from the relevant pump and Stokes 
region, resulting in a poor Raman speetrum. By inereasing the length (L 2 and L 3 bowties) the 
resonanee redshifts and lines up with the pump and Stokes wavelengths. Eor these bowties the 
NTP speetrum starts to emerge. The reported SERS speetra ean henee be attributed to a plasmon 
resonanee effeet sueh that a stable and reprodueible enhaneement faetor ean be assoeiated with 
them, in eontrast to SERS events originating from random surfaee defeets. The inereased overlap 
with the plasmon resonanee, and henee extinetion, also results in a deereased baekground. 

Due to the metal indueed loss, there will exist an optimum number Nopt of patterned antennas 
sueh that the SERS signal is maximized. Sueh an optimum is investigated in Eigurej^for a fixed 
bowtie geometry (a = 60°, L = 100 nm and A = 40 nm) but varying N: N = 10,20,30,40,70 and 
N = 0 whieh is a referenee waveguide. Eaeh waveguide is measured 10 times and the averaged 
Raman speetra are reported in Eigurej^a). The A = 70 signal is not shown sinee it eould not 
be distinguished from the inherent offset signal of the deteetor. Eor a given fixed input power, 
corresponding to roughly 5 mW guided power, the reference waveguide generates a considerable 
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background signal in the 1,340 region, where the strongest NTP peak is expeeted. Funetion- 
alizing the waveguide with inereasing N reduees this unwanted baekground due to the attenuation 
eaused by the nanoantennas. In addition the 1,340 em^^ peak starts to emerge when N inereases. 
The smaller peaks at 1,080, 1,110 and 1,575 em^^ only appear when the baekground is suffieiently 
low. A zoom on the dominant 1,340 em^^ peak (eyan dashed line) is shown in Figure [^b). For 
elarity the baekground is loeally subtraeted. As expeeted, the signal reaehes a maximum value for 
10 < A < 20 and then deeays again with inereasing N. Apart from signal optimization it is however 
equally important to reduee the SiN baekground in order to resolve the smallest speetral features. 
Therefore an analytieal model is developed to outline the interplay between signal enhaneement 
and baekground reduetion, and to derive the relevant figure of merit for on-ehip SERS. 

Figure]^ a) shows a sehematie longitudinal eross-seetion of the ehip eonsisting of N antennas 
spaeed with period A = 10 /im. Eaeh array is eentered on the waveguide with a distanee Li 0.5 
em to the front and baek faeet of the ehip. The NTP monolayer on eaeh antenna will generate a 
forward propagating Stokes power Pa{Xp,Xs) for a given pump power Ppump- This single antenna 
eonversion effieieney Pa{^p,^s) is an antenna dependent faetor ineorporating the integrated field 
enhaneement profile near the metal surfaee and the moleeular density and Raman eross-seetion. 
The total transmission loss indueed by one antenna at wavelength A is given by 1 — whereby 
ex is the linear antenna extinetion. Apart from the intrinsie waveguide losses CUvvgj the pump 
and Stokes light will henee also be attenuated by ep and es respeetively. In the Supplementary 
Information S3 it is then shown that the total Stokes power Ptot generated by an array of N eoated 
antennas is approximately given by 


tot 


- pump 




-2a,,,Li (1-^) 



= FOM{N, Xp, . 


The quantity FOM{N, Xp, Xs) eontains all neeessary parameters to assess the SERS signal strength 
for a given waveguide geometry and is henee eonsidered to be the relevant figure of merit (FOM) 
in eomparing the performanee of integrated antenna arrays. The optimum antenna number Ngpt = 
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log {log (^ 5 ) / log (ep)} / log {esjep). For the particular bowtie antenna studied in Figure the ex¬ 
tinction spectrum e{X) and single antenna conversion efficiency Pa{Xp,Xs) numerically eval¬ 
uated using Lumerical FDTD Solutions (see Methods section). The predicted Nopt for the 1,340 


cm ^ peak is 10 antennas, using the simulated extinctions ep= 1.14 {Ep = 0.58 dB) and es=l .08 


{Es = 0.34 dB), while Pa ~ 2.35 X 10 For each IW of pump power the antenna is therefore 
expected to generate 2.35 fW of guided Stokes power. The Raman enhancement factor is calcu¬ 
lated through EEr = / 3 (Ap)^/ 3 (A 5 )^ in which /3(A) is the local field enhancement (see Methods 
section). In the center of the gap at 5 nm from the tip of the antenna (marked by the black dot in 


Figure [^a)) an EEr 1.42 x 10^ is expected . Apart from the relevant NTP signal, the SiN itself 


generates a considerable background while the pump beam is propagating along the waveguide. 
This background signal P^g can be approximated by 



in which Pb{^p, A 5 ) is a waveguide dependent factor incorporating the specific modal field profile 
and the SiN molecular density and cross-section (see Supplementary Information S4). 

Our analytical model and the associated numerical calculations will now be benchmarked 
against the spectra from Figure to verify whether the theoretically estimated power values corre¬ 
spond to the experimentally obtained absolute Raman powers. To this end, the NTP signal strength 
at 1,340 cm \ obtained from Figure |^b), is analyzed as a function of the antenna number N and 
compared with the theoretical estimations. Furthermore, the background associated shot noise 
is also calculated. The results are depicted in Figure Qb)- While the ideal model assumes N 
identical antennas, the fabricated antennas will show differences among each other resulting in 
changes of ep, es and Pa from one antenna to the other. A generalized model incorporating po¬ 
tential differences in ep, es and Pa is described in the Supplementary Information S5. In order 
to estimate the uncertainty on these experimental parameters a randomized fit to the generalized 
model has been applied. A set of normally distributed numbers is generated for each of the three 
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parameters and then plugged into the generalized model to ealeulate the distribution of signal and 
shot noise eounts, defining an area within whieh the probable signal (blue area) and noise (red 
area) eounts are situated. The mean values of these distributions are extraeted from an initial eon- 
strained fit to the ideal model (dotted lines), while its standard deviations are ehosen sueh that the 
eorresponding signal and shot noise distributions eover all experimental datapoints (red and blue 
dots). Based on the randomized fit it is possible to estimate the spread on the experimental param¬ 
eters: Ep ^ 0.49 ±0.11 dB, Es ^ 0.35 ±0.11 dB and Pa = (2.60 ±0.77) x lO^^^ (theoretieally 
Ep = 0.58 dB, Es = 0.34 dB and Pa = 2.35 x 10“^^ were predieted). The theoretieally predieted 
parameters are all within the error bars of the experimentally fitted data, whieh elearly establishes 
the validity of our model and its ability to provide quantitative predietions of the absolute Raman 
power eoupled into a single mode waveguide. Given this exeellent eorrespondenee, we expeet the 
fabrieated struetures to have a Raman enhaneement faetor EEp on the order of 1.42 x 10"^ near 
the two antenna gap tips. Deereasing the gap size should boost EEp and Pa by another two or 
three orders of magnitude. From the fitting values the optimum antenna number is estimated to 
be 11 ± 3 (eompared to 10 theoretieally). The single antenna eonversion effieieney Pa shows that 
the fabrieated antennas produee (2.60 ±0.77) fW of guided Stokes power for eaeh IW of guided 
pump power. Compared to the free spaee Raman seattering Pq of a single NTP moleeule in a bulk 
air environment Pa ~ 3.94 x lO^Fo- This ineludes the exeitation and emission enhaneement of all 
moleeules in the monolayer as well as the eoupling effieieney to the guided mode. Sinee only half 
of the Stokes light is earried by the forward propagating mode, the total power eoupled into the 
fundamental TE-mode is therefore 8 x 10^/h- 

Our observations also reveal that a minimum number of antennas Nmin is required to generate a 
deteetable signal (marked by the white square in Figure j^b)). If < Nmin then the shot noise still 
dominates on the signal. It has to be noted however that the relevant signal is generated in a very 
small region (A^ — 1)A eompared to the overall length 2Li + {N — l)ARi 2Li, while the shot noise 
is mainly attributed to this non-useful length 2Li. Chip designs whieh allow a separation of the 
signal from the baekground are expeeted to have Nmin = 1 sueh that signals originating from one 
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single antenna ean still be deteeted. As a result it would beeome possible to simultaneously probe 
large areas of analytes (> and deteeting all SERS events, originating from different loeations, 
by monitoring just a single waveguide output in eontrast to mieroseopy based systems where one 
has to serially sean all hotspot loeations. 

The work presented here paves the way towards the effieient design of evaneseently eoupled 
nanoantennas for on-ehip exeitation and emission enhaneement in the 700-1000 nm region. Due to 
the low fluoreseenee, negligible water absorption and the availability of high quality and low-eost 
sourees and deteetors this region is of partieular interest for Raman sensing.l^In eombination with 
other on-ehip speetral funetionalities, sueh as arrayed waveguide gratings, ^the presented platform 
is foreeasted to allow multiplexed deteetion of extremely weak Raman signals on a highly dense 
integrated platform. We also envisage that integrated nanoantennas, similar to the ones reported 
here, eould be used as transdueer between quantum dot emitters and the fundamental waveguide 
mode, potentially enabling applieations in on-ehip quantum eommunieation and quantum eompu- 

tation.l^^^Q] 


Methods 

Fabrication details 

The fabrieation eonsists of a 2-step e-beam lithography proeess. In the first step the nanoplasmonie 
antennas are patterned on top of a slab Si/Si 02 /SiN wafer using a positive PMMA e-beam resist. 
After PMMA exposure, the samples are developed in a 1:1 MIBK:IPA solution after whieh a 2 
nm Ti adhesion layer and 30 nm Au layer are deposited in a eommereial Pfeiffer Spider sputter 
system. The samples are then immersed in aeetone for lift-off. In the seeond step the waveguides 
are defined. After metal lift-off a negative ma-N 2403 resist is spun, exposed and developed in ma- 
D 525. An e-spaeer is also spun on top of the ma-N 2403 to avoid eharging effeets. The developed 
samples are then etehed with an ICP plasma (C 4 F 8 /SF 6 mixture) in a eommereial Oxford Plasmalab 
system. After resist strip and eleaning, the samples are immersed overnight in a 1 mM NTP:EtOH 
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solution and subsequently rinsed with pure ethanol to remove the residual NTP. A self-assembled 
monolayer of NTP is assumed to form on the Au surface through a Au-S bond.^^ 


Numerical Simulations 


Numerical simulations were performed with Lumerical FDTD Solutions. We used a refractive 
index of Urib = 1.9 for the SiN rib (with width w^b = 700 nm and height h^ib = 220 nm), riudad = 
1.45 for the Si02 undercladding and ritdad = 1 for the top cladding (air). The Si substrate was not 
taken into account since the real oxide cladding is thick enough to avoid substantial power leakage 
to the Si such that the numerical results faithfully represent the actual experimental conditions. A 
thin native oxide layer (t„ox = 2 nm) between the SiN and the Ti has also been incorporated.^ 
The metal stack thicknesses are fixed to tn = 2 nm and t^u = 30 nm and a built-in refractive index 
model for Au (Johnson and Christy!^ and Ti (CRC^^ is used. An additional surface layer with 
thickness t^xp = 1 nm and index hmtp = 3 is used to model the NTP monolayer. The antenna 
region (including the Ti adhesion layer and the NTP monolayer) is meshed with a uniform mesh 
of 0.5 nm in the plane of the antenna (yz-plane) and 2 nm in the x-direction. A mesh refinement 
to 1 nm is applied in regions where the thicknesses in the x-direction are < 2 nm. The estimated 
surface area of an NTP molecule is 0.18 nm^, so the surface density is then = 5.56 x 10^® 
molecules/m^.l^ The Raman cross section is (7 ~ 0.358 x 10^^® cm^/sr, which was obtained by 
applying the scaling to the original data of the 1,340 cm^^ line.l^ Single antenna extinction 
spectra £’(A) (in dB) are calculated through E{X) = Trgf{X) — Tant{^) in which Tref{X) is the 
power transmission (in dB) through the reference waveguide and Ta„t{X) the power transmission 
(in dB) of a waveguide functionalized with one antenna. Linear extinction spectra e(A) = ex are 
then given by e{X) = 10^A field and index profile monitor are used to extract the local field 
|E(r, A)| and index n{r) around the antenna. The single antenna conversion efficiency 


Pa{Xp,Xs) 


Z'-VJ-P (^fJ„(r)2\E-(r,Xp)\^dr) (//„(r)2|E'"(rXs)|"<ir) 


8 



is calculated by integrating the local fields over the effective monolayer volume Vm in which the 
index satisfies n{r) = n^^TP (see Supplementary Information S3). The group index of the 
waveguide mode is ng{X). The denominator is calculated using the modal fields E'”(r,A) of a 
non-functionalized reference waveguide and the local field enhancement is given by the ratio of the 
local and modal electric fields: /3(r, A) = At a certain position, the Raman enhancement 

factor EFr is calculated as EFr = j8(Ap)^/3(A5)^. Numerically calculated values are mentioned in 
the main text. 
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Figure 1: Evanescent excitation and collection of SERS spectra, (a) 

Schematic of the chip consisting of single mode SiN waveguides (blue) on an 
Si 02 undercladding (gray), functionalized with an array of gold bowtie anten¬ 
nas (yellow). All antennas are coated with an NTP monolayer (purple dots), 
evanescently excited by the fundamental TE-mode (red). The NTP Stokes sig¬ 
nal (green) is collected by the same mode, (b) Scanning electron microscope 
image of a functionalized waveguide. The white arrows indicate antenna posi¬ 
tions. The inset shows a zoom of a typical antenna, (c) Raman spectra of a 
waveguide functionalized with 40 antennas. The cyan shaded areas mark the 
NTP Stokes peaks while the black dashed lines represent peaks attributed to 
interference effects of the plasmonic array. Before coating, the waveguide al¬ 
ready generates a Raman background in itself (top). After coating, NTP peaks 
emerge (bottom). 


Figure 1 
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Figure 2: Signal dependence on the plasmon resonance.(a) Bowtie an¬ 
tenna geometrical paramters: length L, gap A and apex angle a. (b) Single 
antenna extinction spectra for 3 different bowtie antennas with fixed apex an¬ 
gle a = 60° and gap A = 40 nm but varying length: Li = 90 nm (orange), 
^2 = 115 nm (green), L 3 = 140 nm (blue). The red and cyan shaded lines cor¬ 
respond to the pump and Stokes wavelengths respectively, (c) Corresponding 
Raman spectra of the waveguides functionalized with these 3 bowtie anntenas 
(orange,green,blue) and Raman spectrum of the reference waveguide (red). 


Figure 2 
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Figme 3: Signal dependence on the number of antennas N. (a) Raman 
spectra of a reference waveguide (REF) and waveguides functionalized with 
N = 10,20,30,40 antennas, (b) Zoom on the 1,340 cm“'^ signal peak. The 
background is locally subtracted to obtain the pure NTP signal. 


Figure 3 
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Figure 4: On-chip SERS model and fit to the experimental data, (a) 

Longitudinal cross-section of the chip. The pump beam, with power Ppump at 
wavelength Ap, excites a Stokes signal Ag of which the total power at the output 
facet is Ptot- The plasmonic array consists of N antennas with period A and 
is separated Li from both the input and output facet of the chip. Each of the 
antennas generates Pa guided Stokes power for a given pump power. Apart from 
the waveguide losses a^g, the pump and Stokes light is attenuated by the pump 
ep and Stokes eg extinction respectively, (b) Signal (blue circles) and shot noise 
(red circles) at the 1,340 cm“'^ peak. The dotted lines represent a constrained fit 
to the ideal model while the shaded areas represent a randomized fit to estimate 
the uncertainty on the experimental parameters. The white square denotes the 
minimum number of antennas N^in required to generate a detectable signal. 


Figure 4 
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SI: SERS measurement setup 

SERS spectra are measured with the setup depicted in Figure S1. A tunable Ti:saph laser is set to a 
pump wavelength of Ap = 785 nm (red) after which the polarized beam passes through a half-wave 
plate (A/2) in order to rotate the polarization to a TE-polarized beam. A beamsplitter BSl then 
splits the beam into two parts (solid and dashed red line). The solid path is used to generate the 
forward propagating Raman beam and passes through a laser line filter (FEE) at Ap for side-band 
suppression before it is coupled into the chip by an aspheric lens (ASPH). The output beam is then 
collected with an objective (OBJ) and passes through a polarizer P (set to TE) before it is filtered 
by a dichroic mirror which reflects the pump beam and transmits all Stokes wavelengths (green). 
The Stokes light is collected into a fiber using a parabolic mirror collimator (PMC) after which 
the fiber is split by a fiber splitter (FS) of which 1% goes to a power meter (PM) and 99% to a 
commercial spectrometer from ANDOR (Shamrock 303i spectrometer and iDus 416 cooled CCD 
camera). Mirror M2 blocks the second (dashed) path during the measurement but can be removed 
for alignment purposes. The camera (CAM) and the 1% fiber tap are used during alignment and 
to measure the transmitted power Pj. In order to align the sample we initially set the wavelength 
of the Ti:saph to Xj = 800 nm (such that it can be transmitted through DM1 and collected in the 
power meter) and maximize the transmitted power Pt in the forward path. After optimizing the 
transmission the laser is tuned back to Ap = 785 nm. 
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Supplementary Figure SI: Measurement setup. Ti:saph: tunable Ti:saphire laser emitting the 
pump beam at Xp = 785 nm, CAM: camera, PM: power meter, SR 303i and iDus 416: spectrometer 
and cooled CCD detector from ANDOR, BSl: beamsplitter, A/2: half-wave plate, LLF: laser line 
filter for 785 nm, P: polarizer, M1/M3/M4: fixed mirrors, M2: removable mirror/beam block, 
OBJ: objective (50X, NA=0.9), ASPH: aspheric lens (NA=0.5), S: sample stage, DM1: dichroic 
mirror (reflection R and transmission T shown), PMC: parabolic mirror collimator (EFL=15 mm, 
NA=0.2), FS: fiber splitter. 
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S2: Extinction measurement 


Single antenna extinction spectra, resulting from the plasmon resonance, are measured with the 
setup depicted in Figure S2. Light from an NKT EXR-4 supercontinuum source (SC) is filtered 
through a near-IR filter (NIRF) to filter out the relevant wavelength region. Subsequently it is 
coupled in a fiber with a fiber coupling unit (FC). This fiber is plugged into a fiberbench consisting 
of 3 parts (fixed to the same bench): an achromatic fiber collimator (C) which converts the fiberized 
light to a free-space beam, a free-space broadband polarizer (P) which polarizes the unpolarized 
light into a TE-beam and an aspheric lens (ASPH) used to focus the free-space beam on the input 
facet of the chip. This fiberbench (marked by the light-blue area) is mounted on a piezo-controlled 
stage (XYZ) in order to precisely couple the supercontinuum light into the chip. At the output 
facet a lensed fiber (LF) is used to capture the transmitted light. This lensed fiber is also connected 
to a piezo-controller for accurate positioning. Finally the light is coupled to an optical spectrum 
analyzer (OSA) and the spectra are read out by Python controlled software (PC). 
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Supplementary Figure S2: Resonance measurement. SC: supercontinuum source, NIRF: near-IR 
filter, FC: fiber coupling unit, C: achromatic fiber collimator, P: free-space broadband polarizer, 
ASPH: aspheric lens (NA=0.68), LF: lensed fiber, XYZ: piezo controller stage, S: sample stage, 
OSA: Optical Spectrum Analyzer, PC: OSA control using Python based measurement framework. 
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The single antenna extinction curves E{X) (in dB) can then be calculated through E{X) = {Tref 
TNant)/N in which Tref is the power transmission (in dB) through the reference waveguide and 
Tp^ant the power transmission (in dB) of a waveguide functionalized with N antennas.!^ In Figure 
S3 the single antenna extinction curve of a waveguide functionalized with N = AO bowtie antennas 
(a ~ 60°, // 100 nm and A 40 nm) is shown. The extinction exhibits periodic fringes on the 

broad envelope which are attributed to interference effects of the plasmonic array. On one hand, 
the array forms a multiple Fabry-Perot resonator of which the expected free spectral range FSR of 
19.9 nm around 1,250 cm ^ (using the fabricated array period of A = 10/im) matches well with 
the experimentally obtained value (~ 20.8 nm). On the other hand, the waveguide mode interferes 
with the radiative decay of the plasmon mode. This will affect the specific lineshape and strength 
of the fringes. The spectral positions at which we see sudden changes in the Raman background 
(see main text) coincide with the fringes observed on the extinction curves. Hence these features 
are not attributed to specific Raman lines. 



Supplementary Figure S3: Single antenna extinction spectrum. Single antenna extinction spec¬ 
trum of the waveguide functionalized with = 40 bowtie antennas {a 60°, H ^ 100 nm and 
A 40 nm). 


20 
















S3: Derivation of the analytical on-chip SERS model 

The power /’vvg(ro, coupled into the forward propagating waveguide mode E'”(r, A) as a result 
of a radiating dipole at position Tq is given by 

fivg(rQ,A) _ 3 ng{X) £o£(ro)|e^-E'’^(ro,A)|^ 

fi, «„ \n) /X,^^e„<;(r)|E"(r.A)p*|^.,„ 

where Po — is the power radiated in free space at a wavelength A. — e(r) the relative 

permittivity at position r, ng{X) the group index of the waveguide mode, n,„ the refractive index of 
the medium in which the dipole is placed and the unit vector along the dipole moment vector 
doThe integral in the denominator is calculated over a waveguide cross-section A^g in the xy- 
plane and evaluated at the dipole position z — zo (coordinate axes are defined in Figure S4(a) and 
S4(b)). 

Functionalizing the waveguide with an antenna however introduces a perturbation to the modal 
field E'”(r, A) of the waveguide. In order to investigate whether formula (1) is still valid when 
the waveguide is functionalized with a metallic antenna, two sets of simulations were performed 
. In the first set (Figure S4(a)) the coupling of dipole radiation into the fundamental TE-mode 
is investigated. This allows an explicit evaluation of the left-hand side of equation (1). Dipole 
sources with fixed dipole moment vector do were placed at different positions ro around the antenna 
(shown as black dots in Figure S4(a)). In the second set (Figure S4(b)) the fundamental TE-mode 
is launched into the waveguide functionalized with the same antenna in order to calculate the fields 
at positions rg where the dipole sources in the first simulation set were located. In this way it 
is possible to calculate the right-hand side of equation (1). The integral in the denominator is 
approximated in all calculations by evaluating it on a reference waveguide without antenna, so 
using the true modal fields and not the perturbed fields. This approximation allows a sufficient 
accurate evaluation of the right-hand side of equation (1) as is confirmed in Eigure S4(c) where 
the simulation results are depicted. The red curve represents an incoherent superposition of the 
power coupled into the TE-mode as a result of the complete set of dipole emitters. The dashed 
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blue curve is the predicted power that couples into the TE-mode, and is calculated as an incoherent 
superposition of the predicted coupled powers for each dipole. It is clear that the predicted power 
coupled into the TE-mode (simulation set 2) matches very well with the explicit calculation using 
dipole emitters (simulation set 1). 
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Supplementary Eigure S4: Dipole radiation, (a) Simulation set 1: collection of incoherent dipoles 
emitting into the fundamental TE-mode (red), (b) Simulation set 2: Calculation of the electric 
fields at the dipole positions of simulation set 1 using a fundamental TE-mode excitation. Based 
on these fields one can extract the predicted dipole emission (dashed blue), (c) Simulation re¬ 
sults: incoherent superposition of the TE-coupled power due to all dipoles (red) and predicted 
TE-coupled power using the electric field values (dashed blue). 


Since the right-hand side of equation (1) faithfully represents the power coupled into the funda¬ 
mental TE-mode, it is used to derive an analytical model predicting the total Stokes power coupled 
into the waveguide as a result of an array of N antennas on top of the waveguide (see Eigure 4(a) 
in the main text). Eirst of all we write the dipole strength |doP as a function of the guided pump 
power Ppump in the waveguide 
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in which is the molecular polarizability and E™(ro, A) = ■ E'”(ro, A) the field strength at the 

dipole position. The integral is evaluated at the input Zinput = 0 of the waveguide .1231 The quantity 


t7(ro, Ap, A5) 


ng{Xp)ng{?is)Xl 

|EJ('-o,W|" 

lEi'CroAs)!" 

(//e(r) 

|E«i 

'(r,Ap)pJr) 

(//£« 

|E'”(r,A 5 )pJr) 


( 2 ) 


solely depends on the modal fields and can be evaluated on a reference waveguide. Both integrals 
in the denominator are evaluated along the same cross-sectional area of the reference waveguide (in 
all simulations rim = !)• The power P^g coupled into the fundamental TE-mode due to a collection 
of incoherently radiating dipoles around antenna n{n= 1... A^) is then approximately given by 


pn 

^wg 


- pump 


po{Xs) 


jI ^ .),ro.Ap.As)j3(ro.X,)2/5(ro.As)2rfro (3) 


in which p is the molecular density, ^(As) = Raman cross section, 123 a(Ap) = ocp 

the waveguide loss at the pump wavelength, e(Ap) = ep the linear antenna extinction at the pump 
wavelength, Vm the volume in which the molecules are situated and /3(ro,Ap/5) = /3p/5(ro) the 
field enhancement factor, at the dipole position tq, at the pump (E) and Stokes (S) wavelength 
respectively: 


/3(ro,A) = 


Er(ro,A)| 

|E^(ro,A)| 


whereby E2”^(ro,A) is the local field around the antenna surface, 
rewrite formula (3) as 


For notational simplicity we 


-2^^ ^ Ea(Ap, 

^pump 


in which 

Ra(Ap,A5) = ^ ^ J J Jv Ap)^^(ro,A5)^Jro 


is a dimensionless antenna dependent factor which incorporates the specific field enhancement 
profile near the antenna surface. Formula (3) takes into account the waveguide loss and the loss 
induced by the antennas in front of the n* antenna as a result of which the actual excitation power 
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decays. The signal (at Stokes wavelength A5) still has to pass N — n antennas and propagate along 
a distance Li-\-{N — n)A such that the power reaching the output of the waveguide is given 

by 


pn 

^wg^out 

Ppump 




= Pa(Ap, A5) ( e-«HAi-A)-«.(Ll+iVA)^^^^iV 


— I 

ep 


Since N antennas contribute incoherently to the signal, the total amount of Stokes light at the output 
of the waveguide is given by 
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For the considered waveguide platform it is reasonable to approximate ap ^ as = O^g such that 
the forward propagating Raman power is given by 
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where L = 1L\ + (N —\)A'^ ILi has been used. This result applies to Stokes light co-propagating 
with the pump beam. For a fixed antenna geometry and Xp and Xg, the optimum number of antennas 
Nopt that should be patterned on a waveguide to maximize the SERS signal is given by 


Nopt — 




For the specific case of a monolayer, one has a molecular surface density ps (number of molecules 
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per m^) rather than a volume density p (number of molecules per m^). In this limiting case the 
single antenna conversion efficiency is theoretically defined by 

Pa{^p,'^s) ^ 2 / Ia 

where the integration is now performed over the surface area covered by the monolayer. Nu¬ 
merically one however always needs to model the presence of such a monolayer by introducing a 
finite thickness As a result an effective monolayer volume Kj ~ x is defined. For any 
numerical evaluation the fields are hence integrated over this volume Vm such that 

III ri(^ro,Xp,Xs)l5{ro,Xpf[i{ro,Xsfdro^tmX [ [ ri(^ro,Xp,Xs)[5{ro,Xpfl5{ro,Xsfdro 
J J Jv„ J Ja„ 

The single antenna conversion efficiency is then given by 

Pa{^P,^s) ^ 2I J J Jv 

From a numerical point of view, an effective volume density pe/f = Ps/tm is introduced which is 
then multiplied with the integral over the effective monolayer volume Vm- 
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In order to correlate the Stokes power at the output facet of the chip and the actual detected 
signal counts C 5 , normalized with the integration time T, we need to take into account the charac¬ 
teristics of all optics in our setup. The Stokes power reaching the detector surface is given by 

Ps ^ yout ^ ^ ^ Y>^t ^ ^ [fOMq 

Ppump) 

= f X Tout X if''Pin)) 

= FOM (/■ X Tout 

in which Y and are the coupling efficiencies in and out of the chip respectively, Pin is the 
power just before the input facet of the chip and Tout is the transmission through all optics between 
the output of the chip and the detector surface. The quantity can be written as 

r= Pt/Tpm 

in which Pj is the transmitted power as measured by the power meter (see Supplementary Infor¬ 
mation SI) and TpM is the optical transmission between the output of the chip and the input of the 
powermeter. This transmission is related to Tout by Tout ~ 100 x Tpm x T^pec in which Tspec is the 
transmission between the input slit of the spectrometer and the detector surface. The factor 100 
stems from the fact that the fiber splitter only transmits 1% of the power to the power meter. So 
eventually Ps can be written as Pg = FOM x Pp x lOOT^pec- By using the spectrometer sensitivity 
X (defined as the number of electrons per count) and quantum efficiency QE (A 5 ) (number of elec¬ 
trons per number of incident photons) Ps can be related to the measured signal counts C 5 per unit 
integration time T by 

(^^\( ^ 

\Xs) \QE{Xs) 

Finally one gets 

^ = FOM X imTspec ^ 
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The transmitted power 100/V is about ImW for the given measurement setup. If we consider 
the 1340 cm^^ line then ^ 1 counts/photon and Tspec ~ 0.44 (using data supplied by the 

manufacturer) such that C* ~ 1.94 x 10^^ counts/sec. C*{A,s) can now be used as a conversion 
factor when the experimental signal counts are fitted to the analytical model (and hence the figure 
of merit). In this way one can rigorously quantify the parameters ep, eg and Pa and compare them 
with the theoretical predictions. 
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S4: Derivation of the background signal 

Similar reasonings as in the Supplementary Information S3 can be applied to derive the total back¬ 
ground signal generated along the waveguide. This background mostly stems from the SiN core 
and is hence generated when the pump light propagates along the total length L of the waveguide. 
For a functionalized waveguide the attenuation of the pump and Stokes light due to the antenna 
array also needs to be incorporated. The dipoles which give rise to the background are mainly 
situated in the core, so for a given core cross-sectional area Acore equation (2) has to be integrated 
over the complete waveguide core. This quantity is defined as 


Vcorei^P: 



77(ro,Ap,As)Jro. 


For a given molecular core density pc and core scattering cross section Oc, the total background 
Pbg propagating in the forward direction is then calculated through 
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After some calculation and using ap ^ as = oCwg and L ^ Li/2 (2Li » {N — 1)A), the total 
forward propagating background is given by 
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S5: Randomized fit to a generalized model 

The model derived in the Supplementary Information S3 and S4 assumes identieal antennas. Due 
to fabrieation errors there will always be differenees among eaeh of the antennas in the array. 
These differenees have an impaet on both the extinetion ep and es as well as on the single antenna 
eonversion effieieney Pa{^P:'^s)- While the experimental data deseribed in this Letter ean be fitted 
well to the ideal model, the fit is not perfeet. In this seetion a randomized fit model that takes 
into aeeount the potential deviations among different antennas is outlined. The power generated 
by antenna n is given by: 

^pump 

Sinee Li >> NA this ean be simplified to 
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The formula however assumes a eonstant pump extinetion from the previous n—l antennas. Fur¬ 
thermore the faetor A 5 ) is also assumed eonstant for eaeh antenna. If we allow that eaeh of 

the antennas has a different extinetion and antenna dependent faetor /’“(Ap, A 5 ) (m = 1 .. .A^), 
then the power generated by antenna n is given by 
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Applying a similar reasoning to the Stokes light that has to propagate along N — n other antennas 
(and a length Li), the Stokes power (generated by antenna n) reaehing the output is given by 
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The total generated Stokes power from all N (potentially different) antennas is then finally 
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For the background one can analogously write 
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(5) 


In this case Pb is a constant since it only depends on the waveguide parameters. The shot noise 
associated to the background is then simply proportional to ^/Pbg■ The number of counts/sec can 
be obtained by multiplying equations (4) and (5) with C*(As) (see Supplementary Information S3). 
Based on the conversion factor C*(A 5 ) and a fit to the experimental data, one can quantitatively 
determine a value for the parameters in the model. For a fixed Ap and Ag, Pa(Ap, As)C* (As) will be 
denoted by P^. 

In order to make a randomized fit of the signal and shot noise at the same time, N random 
numbers (based on a normal distribution) are independently generated for the pump and Stokes 
extinction and for the single antenna conversion efficiency (so three random numbers are generated 
for each antenna in the array). Using these random numbers, formulas (4) and (5) are evaluated. 
This process is repeated 1000 times in order to generate a statistically relevant distribution of the 
possible signal and shot noise counts. Ultimately the mean value /i and standard deviation a from 
the obtained signal and shot noise distributions are extracted for each number of antennas N. The 
Sc-intervals for both distributions ([/i — 3(7, /i + 3(7]) define an area that marks the possible signal 
and noise counts for a given uncertainty on the antenna parameters (according to formulas (4) and 
(5)). These areas are plotted in Figure 4(b) of the main text. 

The mean value of each of the three normal distributions (for ep, es and Pa) is determined by a 
fit of the experimental data to the ideal model. One could obtain a perfect fit to either the signal or 
shot noise data by fitting only one of the two equations ((4) or (5)). The fitting parameters would 
however not generate a perfect fit to the other non-fitted equation since both equations depend on 
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ep and es- Therefore a constrained fit, in which one minimizes the fitting error to both equations 
simuitaneousiy, is performed such that both the signai and shot noise data are represented weii 
using the fitted vaiues for ep and es- The vaiues ep, es and (mind that and not Pa is used 
because we fit the counts/sec) obtained from the constrained fit are then chosen to be the mean 
vaiues (Pep, Pes and pp*) of their respective normai distributions. The standard deviation Ox is 
chosen such that the 30x intervais represent reaiistic deviations from the mean vaiue (x denotes 
one the three parameters). For the extinction this means e.g. that Pepj^ — > 1 (since 1 is 

the iower boundary for the iinear extinction). Since is aiways positive, pp* — 3Op* > 0 shouid 
also be satisfied. For the experimental data shown in Figure 3 of the main text, pgp = 1.1182, 
Peg = 1.0845 and pp* = 5.0394 has been derived from the constrained fit. The sigma value for the 
extinction is chosen to be c = 0.028 such that Ep ^ 0.49 ± 0.11 dB and Eg ~ 0.35 ±0.11 dB. For 
we similarly get P^ ^ 5.04 ± 1.5 such that the single antenna conversion efficiency for the 1340 
cm~^ line is /^ = (2.60 ±0.77) x 10~^^. All of these values satisfy the constraints mentioned 
above and represent realistic deviations of the three fitting parameters. 
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